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i!S ■ ABSTRACT 

Oh' 

O, 

^ . It is shown that during the motion of two initially gravitationally bound 

j^ i spherical shells, consisting of point particles moving along ballistic trajectories, 

one of the shell may be expelled to infinity at subrelativistic speed Vexp < 0.25c. 
'^>( I The problem is solved in Newtonian gravity. Motion of two intersecting shells in 

H i the case when they do not runaway shows a chaotic behaviour. We hope that 

this toy and oversimplyfied model can give nevertheless a qualitative idea on the 
nature of the mechanism of matter outbursts from the dense stellar clusters 

Subject headings: black holes — ejection:active galactic nuclei: stars: chaos. 



1. Introduction 



Dynamical processes around supermassive black holes in quasars, blazars and active 
galactic nuclei (AGN) are characterised by violent phenomena, leading to formation of jets 
and other outbursts. Jet formation is usually connected with the processes in the magnetized 
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accretion disks (Lovelace, 1976; Bisnovatyi-Kogan & Blinnikov, 1976). Formation of qua- 
sispherical outbursts, which are probably observed in quasars with broad absorption lines, 
could be connected with another mechanism. Here we consider a possibility of a shell out- 
burst from a supermassive black holes (SBH) surrounded by a dense massive stellar cluster, 
basing on a pure ballistic interaction of gravitating shells oscillating around SBH. 

Ballistic ejection may be responsible for appearance of stars and observed SN in the 
intergalactic space, as well as the existence of stars between galaxies, suggested by many 
authors. 

Investigation of spherical stellar clusters using shell approximation was started by Henon 
(1964), and than have been successfully applied for investigation of the stability Henon 
(1973), violent relaxation and collapse (Henon 1964; Gott 1975), leading to formation af a 
stationary cluster. Investigation of the evolution of spherical stellar cluster with account of 
different physical processes was done on the base of a shell model in the classical serie of pa- 
pers of L. Spitzer and his coauthors (Spitzer & Hart 1971a,b; Spitzer & Thuan 1972; Spitzer 
& Shapiro 1972; Spitzer & Chevalier 1973; Spitzer & ShuU 1975a,b; Spitzer & Mathiew 
1980). 

Numerical calculations of a collapse of stellar clusters in a shell approximation (Yangu- 
razova & Bisnovatyi-Kogan 1984; Bisnovatyi-Kogan & Yangurazova 1987,1988) had shown, 
that even if all shells are initially gravitationally bound, after a number of intersections 
some shells obtain sufficient energy to become unbound, and to be thrown to the infinity. In 
the Newtonian gravity the remnant is formed as a stationary stellar cluster, and in general 
relativity SBH may be formed as a remnant. 

Formation of outbursts due to the ballistic interactions may happen as a result of in- 
tersection of shells, oscillating around SBH. In the smooth cluster with or without SBH in 
the centre stars evaporate, due to stellar encounters mainly with small kinetic energy, and 
formation of rapidly expelling stars is about 100 times less probable due to predominance of 
gravitational encounters with small momentum transfer Ambartsumian (1938). If the cluster 
is strongly aggregated, consisting of few compact pieces, the encounters between these pieces 
are quite different, collision with large momentum transfer are becoming probable. In this 
case gravitational interaction between compact pieces may lead to the outburst with large 
velocity, and when such event happens in the vicinity of SBH the velocity could become a 
considerable part of the light velocity c. 

Such kind of situation in principle may be produced as a result of galactic collisions 
with close encounter of nuclei, when one of the nuclei strips the matter from the companion 
in the form of a collapsing shells. Interaction of such shell with the stellar cluster may lead 
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not only to collapse into SBH, but also to the reverse phenomena: expelling of the shell with 
a speed, much larger than the average speed of stars in the cluster. The shell is not falling 
into SBH due to high angular momentum of its stars. 

Even more important example of a quasi-spherical mass ejection is a relativistic collapse 
of a spherical stellar system, which is considered (Zeldovich & Podurets, 1965; Lightman & 
Shapiro 1978; Ipser 1980; Rasio et al. 1989) as the main mechanism of a formation of super- 
massive black holes in the galactic centers. Approximation of such collapse by consideration 
of spherical shells is the simplest approach, which reflects all important features of such 
collapse (Gott 1975; Bisnovatyi-Kogan & Yangurazova 1987, 1988). Our consideration is 
related to the motion of stars (shells) which remain outside the newly formed supermassive 
black hole. 

Here we consider a simplifled problem of a motion of two massive spherical shells, each 
consisting of stars with the same speciflc angular momentum and energies, around SBH. The 
aim of this paper is to give an elementary treatment of the ballistic ejection, and to estimate 
the maximal efficiency of this ejection. Note, that in a more complicated case of numerous 
shell intersection this elementary act is a key process of the energy exchange between stars 
and of matter ejection. This is also the elementary process leading to the violent relaxation 
of the cluster Lynden-Bell (1967), studied in the shell approximation by Gott (1975) and 
Yangurazova & Bisnovatyi-Kogan (1984). Chaotization of the motion of two gravitating 
intersecting shells appears as a by-prodiuct of our consideration. The appearence of chaos in 
such a system, which is described by a set of simple algebraic equations, is a rare example. 
For the oversimplified case with a pure radial motion and refiecting inner boundary the 
chaotic shell motion was found by Miller and Youngkins (1997). Here we consider realistic 
shells consisting of stars moving along eliptical trajectories, and their chaotic properties 
would be examined in more details in another work. We find conditions at which one of two 
shells is expelling to infinity taking energy from another shell. We find a maximum of the 
velocity of the outbursting shell as a function of the ratio of its mass m to the mass M of 
SBH using only Newtonian theory of the shell's motion. 

However we introduce cutoff fixing the minimal radius of the potentially outbursting 
shell Tm on the level of few Vg = 2GM /c^. We show, that for equal masses of two shells the 
expelling velocity reach the value Vmax ~ 0.3547fp at m/M = 1.0 , and Vmax > O.Svp was 
obtained at the masses of shells 0.25 -r- 1.5M, where the parabolic velocity of a shell in the 
point of a smallest distance to the black hole Vp may be of a considerable part of c. 

In sections 2,3 we describe outburst effect and in sections 4 we present the evidence of 
the chaos in the system of intersecting shells. The exact solution of these problems in the 
context of General Relativity will be presented in the subsequent paper. 



-4- 



2. Two shells around SBH. 

Physically the nature of the ballistic ejection is based on the following four subsequent 
events. The outer shell is accellerated moving to the center in a strong gravitational field of 
a central body and inner shell. Somewhere near the inner minimal radius of the trajectory 
shells intersect. After that the former outer shell is deccellerating moving from the center 
in the weaker gravitationsl field of only one central body. The second intersection happens 
somwhere far from the center. That may result in the situation when the total energy 
(negative) of the initially gravitationally bound outer shell is becoming positive as a result 
of two subsequent intersections with another shell. The quantitive analysis of this procrss is 
done in this section. 

Equation of motion of a shell with mass m and total conserved energy E in the field of 
a central body with mass M is 

mf ^ Gm{M + m/2) J'^m 
^ = -2 r + ^^' ^^^ 

where v = dr/dt is the radial velocity of the shell and J^m/2r^ is the total kinetic energy 
of tangential motions of all particles, the shell make-up from. The constant J > has 
that interpretation that Jm is the sum of the absolute values of the angular momenta of all 
particles. The term m/2 in (1) is due to the self-gravity of the shell. 

Let's consider two shells with parameters mi, Ji and m2, J2 moving around SBH with 
mass M. Let the shell "1" be initially outer and the shell "2" be the inner one. Then 
equations of motion are: 

"^1^1%) Gm^{M + mi/2 + m2) , Jf^i ,^, 

^i(o) = —^ + ^^, (2) 



^2(0) = — 7, z + ^^ir- \^) 



^ 2^^2(0) _ G'm2(M + m2/2) J|m2 
2 r ^ 2r2 

By the index (0) we mark the initial evolution stage before the first intersection of the shells. 
Assume that both shells are moving to the center and v\,q-. > V2/Q^ . Such shells intersect each 
other at a some radius r = ai and at some time t = ti after wich the shell " 1" becomes inner 
and shell "2" - outer. The equations of motion of the shells during this new stage (designate 
it by the index (1)) are: 

mivl^^) G'mi(M + mi/2) Jfm^ 

Em = —2 + ^^' (4) 

"^2t^2{i) G'm2(M + m2/2 + mi) J|m2 . , 

^2(1) = —^ + ^^- (5) 
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The crucial point now are the matching conditions at the intersection point r = ai, t = ti 
where from one can obtain the initial data to the equations (4), (5) in order to define uniquely 
the evolution during this new stage. In Newtonian theory these conditions are: 

-^1(0) + -^2(0) = -^1(1) + -^2(1); 

^^i(o)(ii) =^i(i)(ii); V2i0)ih) =V2{l)itl), (6) 

i.e. the conservation of the total energy of the system and continuity of the velocities through 
the intersection point. Expressing t>i(o)(ti), f2(o)(^i) from eqs. (2), (3) and t>i(i)(ti), f2(i)(ti) 
from eqs. (4), (5) and equating them following conditions (6) we get: 

_ Grmm2 _ _ Gmim2 
-c^l(l) — -c^l(O) H , -£^2(1) — -£^2(0) • [') 

0,1 (1\ 

At some point r = 02, t = t2 can happen the second intersection after which the shell "1" 
again becomes the outer and shell "2" becomes inner one. It is a simple task to write the 
equation of motion during this third stage (we designate it by the index (2)) and repeat the 
same procedure of matching velocities but now at point r = a2, t = t2. The result is: 

^ ^ Gmim2 ^ „ / 1 1 \ , , 

£^1(2) = ^i(i) = £^1(0) + Gmim2 , (8) 

02 V"^! "^2/ 

^2(2) = ^2(1) H ^^ = ^2(0) - Gmim2 i ) . (9) 

0,2 \fll ^2/ 

Let us describe the situation, when one shell is ejected to infinity after intersection of to 
initially bound shells. We consider a case when 02 is larger then ai so, that the second term 
in (8) has larger absolute value, then the first one, the first shell gains a positive energy 
and goes to infinity. Both shells have initial negative energies -Ei(o) and -^2(0) 5 but with 
small enough absolute values. The first shell takes the energy from the second one, which is 
becoming more bound with larger absolute value of the negative energy -^2(2), according to 



3. Numerical solution. 

Let's illustrate the foregoing scenario by an exact particular example of two shells of 
equal masses. We choose parameters in the following way: 

mi = 7712 = m, E^o) = ^2(0) = 0, Ji < J2, (10) 
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In fact such exact solution represents the first approximation to the more general situation 
when -Ei(o) and -^2(0) are nonzero (negative) but small in that sense that both modulus |-Ei(o)| 
and 1-^2(0)1 S'-re much less than Gvn^ jai. 

We assume that initially both shells are moving towards the center. It follows from 
eqs. (2), (3) that under condition (10) such shells will intersect inescapably at the point 
r = ai, t = ti. After the second intersection at r = 02, t = t2 the shell "1" will be thrown 
to infinity with expelling velocity Viexp- It follows from (8) that 




Viexp= Vi(2)\^^^= \2Gm [—- — ]. (11, 



In order to construct a solution with maximal possible f lexp we consider a case when initially 
inner shell "2" reaches the inner turning point at minimal possible radius r = r^ (see 
introduction), and intersects with the initially outer shell after, during its outward motion. 

To find the second intersection point we need to solve equation of motion (4), (5) for 
shell's trajectories after the first intersection. First, let us consider the bound shell "2". We 
apply eq. (3) to the point r = Vm which gives (see Fig. 1.) 

J2 = G'r„,(2M + m). (12) 

Substituting this expression for J2 into (5) and taking into account (7), (10) we get a 
solution of the equation (5) in the following parametric form: 

riv) =&2(1)(1- 62(1) COST]), 

I ^) (13) 

^^^^ = V g(^ + 3W2) ^'''^^^^""'^'^" 

where the constants 62(1) and 62(1) are: 



M + 3m/2 I 4m(M + m/2)r„ , , 

^'^'^ = 2m "^' "^(^) = V (M + 3m/2)2 V ^ ^ 

on the choice of the constant 6t see below. 

Between the first and second intersection there exists the inner turning point of the 
shell "1" (now inner shell). We take that additional restriction that the shell "1" reaches 
this turning point also at the minimal possible radius r = r^- It is easy to show that a wide 
class of solutions with such restriction really exists. At the point r = r„j equation (4) gives: 

J2 = Gr„X2M + m) + 2Gmrl^a^\ (15) 
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We substitute this expression for Jf into eq. (4) and again with conditions (7) and (10) 
obtain the following parametric form of the solution of this equation: 

rJ O = &i(i)(ei(i)C osh^ - 1), 

I C (16) 

where the constant 6i(i) and ei(i) are: 

M + m/2 2m rm .^„. 

2m M + m/2 ai 

The first intersection point r = ai, t = ti corresponds to some parameters rj = rji, C, = C,i- 
At this point the time t{rii) calculated from (13) should coincide with time t(^i) calculated 
from (16). This fix the constant St'. 



2(1) 

G(M + 3W^(^^-^^«^^"^^^ 
(ei(i)smh^i -^i). 



G(M + m/2) 



Also the values of r{rii) from (13) and r{S,i) from (16) should be the same and equal to ai. 
From this follows parameters rji and ^i: 

coshCi = ^i^i^^^ (19) 

01(1)61(1) 

^2(1) - fll lc.r,^ 

COS rji = —^ (20) 

02(1)62(1) 

The second intersection point r = 02 corresponds to parameters 1] = ri2, C = ^2- The 

application of the analogous coincidence conditions to the equations (13) and (16) at this 
second point gives: 

02 = &i(i) (ei(i) cosh ^2-1) (21) 

02 = &2(i)(l - 62(1) cos 772) (22) 



1(1) 

(6i(i) sinh^ - 6 - ei(i) sinh^i + ^1) 



GjM + m/2) 

2(1) 

-(^2 - 62(1) sin?72 - 771 + 62(1) sin?7i). 



G(M + 3m/2) 



From (14), (17) and (19)-(23) it is easy to see that all these relations determine 02/01 as a 
function of ai/r^ and m/M . 

We introduce now the "parabolic" velocity Vp of the any outer shell at the point of its 
minimal distance to the center r = r„, as: 



V, = !^(^1±M!) , (24) 



Than we have 

7)1 / m I r .. r .. \ 

(25) 



^lea^p / ^^ / f'm ^v 



Vp V M + 3m/2 \ ai 02 

Consequently this ratio is also a function of fli/r^ and m/M. We find numerically the value 
of ai/rjn maximizing the ratio Viexp/vp which value is a function of m/M only. Dependencies 
for Viexp/vp, 02/01, and Oi/r^ as functions of the parameter m/M for these maximizing 
solutions are given in Fig. 2-4. 

The numerical calculations (illustrated by the Fig. 2) shows that the runaway velocity 
Viexp reaches its maximal possible value at m/M k, 1.0 and it is Viexp ~ 0.3547fp. If we 
consider the shells around SBH then the minimal radius Vm of the shell " 1" orbit cannot be 
less then the two gravitational radiuses 2rg = 4G{M + 3m/2)/c^. In the extreme case when 
Tm = 2rg we have Vp ~ c/\f2 and for the maximal possible runaway velocity in the family of 
solutions, characterized by (12), (15) we get Vx^xp ~ 0.25c. 



4. Chaos in the shell motion. 

The first evidence that the motion of two intersecting shells can show chaotic character 
was given by B.N. Miller & V.P. Youngkins (1997). They investigated the special case 
when the central body is absent and particles consisting the shells are moving only in radial 
direction (in our notation M = and Ji = J2 = 0). This situation, however, cannot model 
astrophysical cluster with massive nuclei, and also the problem of the influence of central 
Newtonian singularity arise which need some additional care. In any case a study of more 
physically realistic models with nonzero M, Ji and J2 from the point of view of possible 
chaotic behaviour represents essential interest. We report here some results for such more 
general two shells model which was investigated in the previous sections but again for the 
shells with equal masses. We consider now only oscillatory regime of motion without any 
runaway effects. 

The shell motion in the Newton gravitational field is completely regular, but at presence 
of intersections the picture changes qualitatively. The shell intersections result to chaos in 



-9- 



their motions. Character of chaos depends, mainly, on mass ratio of a shell and a central 
body. At a small shell mass the motion of the shells occurs basically in the field of the 
central body, and after one intersection there is a little change in a trajectory of each shell. 
Consider a case when angular momentum parameters J and energies of shells are close 
to each other, and as a result of two intersections energy is transferred to the shell "1". 
Since shell parameters changes are not large, the mutual positions and speeds of shells 
after intersections do not change strongly. Then at the next intersection the part of energy 
will again be transferred to the shell " 1" . After a large number of intersections, however, 
the mutual positions changes essentially. Then the energy starts to flow to the shell "2", 
consequently we observe beating in oscillation of shells. In the case when the masses of the 
intersecting shells are large enough, the energy is transferred very intensively between them. 
It changes trajectories of shells essentially so that the next intersection does not resemble at 
all the previous one. In this case a shell's motion at once gets properties of a randomness. 

Shell oscillations with mass ratio m/M = 0.015 are presented in Fig. 5. In the case of 
massive shells the exchange of energy between shells occurs more intensively. As a result 
we have the obviously expressed chaotic behaviour of shells. Chaotic shell oscillations with 
m/M = 0.08 are presented in Fig. 6. In Fig.6 a full randomness of behaviour of the shells is 
seen obviously. 

Let us note, that when the exchange of energy between shells is sufficiently large the 
strong randomness appears. Shells may exchange their positions during the motion. The 
external shell may become inner one and vice versa. One of the shells may give up practically 
all its energy of the orbital motion, which it may loose at constant angular momentum, i.e. 
the orbits of its particles become almost circular. It may be seen in Figs. 6, 7 at those moments 
when the orbit of the inner shell fills a narrow band by itself. 

Let us note that in our calculations the motion of the shells is very sensitive to initial 
parameters, what gives an additional evidence that intersecting shells represent a really 
chaotic system. When we change only precision of integrator from 10~® to 3 x 10~^, after 
several intersection motion of shells becomes absolutely different from its initial behaviour in 
Fig.6 at the same initial conditions. If we change the angular momentum parameter of one 
of shell from 1.40 to 1.41, then in the same way as in previous case, motion of shells becomes 
absolutely different after several intersection. We can make statement that the behaviour of 
shells is unstable and is strongly affected by any perturbation in initial parameters of shell 
motion. 

The example of a chaotic behaviour of two intersecting selfgraviting shells, moving in 
their own gravitational field, without a central mass, is represented in Fig. 7. 
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Chaotic behaviour of the intersecting selfgravitating shells in the Newton gravity is a rare 
example of a chaotic dynamical system, described fully analytically by algebraic relations. 



5. Discussion 

We have shown that pure gravitational interaction of two spherical shells and central 
object may lead to considerable energy exchange and ejection of one of the shell with a 
subrelativistic speed ~ 0.25c. Motion of stars around SBH may lead to similar effects due 
to their interaction, and we may expect formation of very rapidly moving stars, which could 
overcome the galactic gravity and runaway into the intergalactic space. In this way very 
extended galactic halo may be formed. Ballistic mechanism of energy exchange between 
gravitating particles may be important in the evolution of spatial distribution of cold dark 
matter, when gravitational instabilities are developing in masses, much less than galactic 
ones, and gravitationally bound massive objects, consisting mainly from the dark matter 
interacting only due to its gravitation, may be formed. Interesting examples of ballistic 
interaction of satellites with planets and stars, their acceleration and ballistic control are 
presented in the book of V.V. Beletskyi (1977). 
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Figure 1. Time dependence of radiuses r of two shells (in units r^) on time t (in units 
r„Jvp, Vp = a/2G'(M + 'im/2)/r„i)- The shells intersect at points ai = 1.1401, 02 = 3.6316. 
Here m/M = 0.1, and after the second intersection the first shell is running away with 
velocity at infinity Viexp = 0.23415fp. 

Figure 2. Dependence of the maximum runaway velocity viexp in units Vp on the ratio 
of the shell mass to the mass of central body m/M. The maximum value of Viexp = 0.3547 
corresponds to m/M = 1.0. 

Figure 3. Dependence of the ratio of intersection radiuses 02/01 on m/M for the solu- 
tions with maximum runaway velocities. 

Figure 4. The value ai/r„j (the ratio of the radius of the first intersection ai to the 
minimal radius of the approach to the central body r^) as a function of m/M for solutions 
with maximum runaway velocities. 

Figure 5. Shell oscillations with mass ratio m/M = 0.015. 

Figure 6. The chaotic motion ri 2(t) of shells with masses m/M = 0.08. 

Figure 7. The chaotic motion ri^2{t) of two selfgravitating shells in their own gravita- 
tional field without a central mass. 
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Fig. 1. — Time dependence of radiuses r of two shells (in units r™) on time t (in units 
'Tm/vp., Vp = a/2G'(M + 3m/2)/rm). The shells intersect at points ai = 1.1401, 02 = 3.6316. 
Here m/M = 0.1, and after the second intersection the first shell is running away with 
velocity at infinity Viexp = 0.23415fp. 
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Fig. 2. — Dependence of the maximum runaway velocity fiexp in units Vp on the ratio of 
the shell mass to the mass of central body m/M. The maximum value of Viexp = 0.3547 
corresponds to m/M = 1.0. 
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Fig. 3. — Dependence of the ratio of intersection radiuses 02/01 on m/M for the solutions 
with maximum runaway velocities. 
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Fig. 4. — The value oi/r^ (the ratio of the radius of the first intersection ai to the minimal 
radius of the approach to the central body r„j) as a function of m/M for solutions with 
maximum runaway velocities. 
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Fig. 5. — Shell oscillations with mass ratio m/M = 0.015. 
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Fig. 6. — Chaotic shell oscillations with mass ratio m/M = 0.08. 



-19- 



120 



100 




xlO 



Fig. 7. — The chaotic motion ri 2(t) of two selfgravitating shells in their own gravitational 
field without a central mass. 



